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Abstract 

In this paper, a supersymmetric extension of the polytropic gas dynamics equations is 
constructed through the use of a superspace involving two independent fermionic variables 
and two bosonic superficlds. A superalgebra of symmetries of the proposed extended model 
is determined and a systematic classification of the one-dimensional subalgebras of this su- 
peralgebra is performed. Through the use of the symmetry reduction method, a number 
of invariant solutions of the supersymmetric polytropic gas dynamics equations are found. 
Several types of solutions are obtained, including algebraic-type solutions and propagation 
waves (simple and double waves). Many of the obtained solutions involve arbitrary functions 
of one or two bosonic or fermionic variables. In the case where the arbitrary functions in- 
volve only the independent fermionic variables, the solutions are expressed in terms of Taylor 
expansions. 
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I Introduction 



Symmetries of fluid dynamics have recently attracted a considerable amount of interest 
among physicists and mathematicians. Dispersionless hydrodynamic systems can be stud- 
ied from a variety of points of view, of which the best understood are (1 + l)-dimensional 
systems. Such systems include, among others, polytropic gases, the Chaplygin gas, the 
Born-Infeld equations and hydrodynamic systems written in terms of Riemann invariants. 
Such systems also include strings and Nambu-Goto membranes (see e.g. R. Jackiw^ and 
references therein). 

The purpose of this paper is to construct a supersymmetric version of the polytropic 
gas dynamics equations in (1 + 1) dimensions 

Pt + PxW + pw x = 

(1) 

w t + ww x + Ap 1 p x = 0, 

where w is the velocity of the fluid, p is its density and 7 is the polytropic exponent. 
These equations are linked with the ones given in Das and Popowicz^ (p. 3, equation (1)) 
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dt K )XJ dt V 2 (7-1)/* 

(given in the form of conservation laws) through the transformation w = —u, p = v. 
For the case where the polytropic exponent is 7 = —1, we get the Chaplygin gas. Su- 
persymmetric versions of the Chaplygin gas in (1 + 1) and (2 + 1) dimensions were for- 
mulated through the parametrization of the action for a superstring and supermembrane 
respectively.^"^ The case 7 = 2 refers to the dispersionless limit of two bosonic equa- 
tions.^ When 7 = 3, the equations represent the dispersionless limit of two non-interacting 
Korteweg-de Vries equations (known as Riemann invariants equations) DSE Finally, the 
7 = 4 case is the dispersionless limit of a Boussinesq equation!^ 

In the context of hydrodynamic-type equations, other models have been supersym- 
metrized and analyzed from the group-theoretical point of view. Examples of such models 
include the Korteweg-de Vries equationp^J] the Kadomtsev-Petviashvili equation,^ the 
scalar Born-Infeld modelp^ a Gaussian fluid flow^l and a hydrodynamic system expressed 
in terms of Riemann invariantsP^"^ The results are of interest since we can construct 
certain classes of solutions with the freedom of arbitrary functions of one or two arguments 
involving bosonic and fermionic variables. 

In the case of several integrable systems, the Lax pairs and multisolitonic solutions 
have been found for their supersymmetric extensions (see e.g.^^). It is natural in this 
context to consider other types of wave superpositions. It has been established that 
classical hydrodynamic-type systems admit simple and multiple Riemann wave (scattering 
and non-scattering) solutionsP^ The question arises as to whether or not the same holds 
true for a supersymmetric extension of such a model. We will answer this question in the 
case of the polytropic gas dynamics equations (JT|) by performing a symmetry analysis. 

The supersymmetric extension that we construct in this paper differs from that of 
Das and Popowicz^ in the sense that our superspace includes two independent fermionic 
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variables instead of one. Also, the superfields which we propose are bosonic-valued instead 
of fermionic- valued. This is a different form of supersymmetric extension for the equations 
(PQ) which, to our knowledge, has not previously been considered. We perform a systematic 
group-theoretical analysis of our constructed supersymmetric model. 

This paper is organized as follows. In section 2, we construct a supersymmetric ex- 
tension of the polytropic gas dynamics equations (jTJ through a superspace and superfield 
formalism involving two independent fermionic variables. In section 3, we present in detail 
a number of symmetries of our supersymmetric system. For this purpose, we use a general- 
ization of the method of prolongation of vector fields, extended so as to include both even 
and odd Grassmannian variables. The symmetry criterion is adapted to this case, and 
the superalgebra associated with the supersymmetric system is determined. We classify 
the one- dimensional subalgebras of this superalgebra into conjugation classes under the 
action of the associated supergroup, and this allows us to perform symmetry reductions 
in a systematic way. In section 4, we proceed to find the invariants and reduced sys- 
tems corresponding to those subalgebras with standard invariants structures. The other 
subalgebras, possessing a nonstandard invariant structure, are discussed separately. We 
propose a new method for solving the reduced equation. This leads to solutions of the 
supersymmetric polytropic gas model, many of which contain the freedom of arbitrary 
functions of one or two arguments which are expressed in terms of bosonic or fermionic 
variables. In the case when the arbitrary functions involve only the independent fermionic 
variables, we were able to express the solution in terms of a naturally truncated Taylor 
expansion. Finally, section 5 contains the final remarks and future outlook. 



II Supersymmetric extension 

The supersymmetric extension is constructed by first considering the following superspace 
and superfield formalism. We extend the space of independent variables {(x,t)}, which 
presently includes only bosonic variables, to the superspace {(x, t, 6\, 6*2)} which also in- 
cludes the independent fermionic variables 9\ and 9 2 . The variables x and t represent 
the bosonic (even Grassmannian) coordinates on 2-dimensional Minkowski space, while 
the quantities 9\ and 9 2 are anticommuting fermionic (odd Grassmannian) variables. We 
replace the real bosonic-valued fields of velocity w(x, t) and density p(x, t) by the bosonic 
superfields defined as 

W(x,t,e 1 ,9 2 ) = w(x,t) + 9 1 4> 1 (x,t) + 9 2 cf> 2 (x,t) + 9 1 e 2 F(x,t), 

(3) 

P (x, t, 0i, 2 ) = p(x, t) + 6xfa(x, t) + 9 2 ip 2 (x, t) + 9 1 9 2 G(x, t) 

respectively, where 0i, <fi 2 , ipi and ip2 are fermionic fields and F and G are bosonic fields. 
The supersymmetric extension of the polytropic gas equations JT]) is constructed in such 
a way that it is invariant under the supersymmetry transformations 

x-tx-rijOi, 9 1 -*8 1 + ri 1 and t -> t - 77^2, $2^02+%, (4) 

where 771 and rj 2 are fermionic constants. In what follows, we use the convention that un- 
derlined constants represent fermionic parameters. These transformations are generated 
by the infinitesimal supersymmetry generators 

Qi = 9 01 - 8id x and Q 2 = d 02 - 9 2 d u (5) 
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which satisfy the anticommutation relations 



{QuQi} = -20*, {Q 2 , Q 2 } = -2dt. (6) 

Here, we use the concept of supercommutation in the sense that we take the commutator 
of either two bosonic quantities or of a bosonic quantity and a fermionic quantity, but 
the anticommutator of two fermionic quantities. In order to make our superfield theory 
manifestly invariant under the action of the supersymmetry generators Q x and Q 2 , we 
write the supersymmetric system in terms of the covariant derivative operators 

D 1 = d ei +9 1 d x and D 2 = de 2 + 9 2 d tl (7) 

which possess the property that they anticommute with the supersymmetry generators 
Qx and Q 2 

{D x ,D x }=2d x , {D 2 ,D 2 } = 2d t , 
{£>!, D 2 } = {D x , Q x } ={D 1 , Q 2 } = {D 2 , Q x } = {D 2 , Q 2 } = 0. U 

We supersymmetrize the polytropic gas dynamics equations ([I]) directly at the equation 
of motion level instead of from its Lagrangian form. In order to build the supersymmetric 
extension, we evaluate the covariant derivatives of the superfields W and P of various 
orders. The most general supersymmetric extension of the polytropic gas equation (pQ) is 
constructed by considering linear combinations of the products of the various covariant 
derivatives of the superfields W and P. These multiply together to produce the given 
terms as coefficients whose component reproduces each term of the classical equations 
(JTJ). The result of this analysis gives the following form of the supersymmetric polytropic 
gas dynamics equations 

P t + a(D x D 2 P) x W + b{D 2 P) x {D x W) + c(D x P) x (D 2 W) 

+ ( c _ b _ a _ l)P x (D x D 2 W) + dP{D x D 2 W) x + e{D x P){D 2 W) x (9) 

+ f(D 2 P)(D x W) x + {e-f-d- 1){D X D 2 P)W X = 

and 

W t + gW(D x D 2 W) x + h(D x W)(D 2 W) x + t(D 2 W)(D x W) x 

+ (h-i- g - 1){D X D 2 W)W X + j{D x D 2 PyP x + kiD.D.Py-^D^iD.P), 

+ l{D 1 D 2 Py-\D l P){D 2 P) x + m{D x D 2 Py- 1 P{D l D 2 P) x ( } 

+ ((-iy +1 A + J + k + m-l) (D 1 D 2 Py- 2 (D 1 P)(D 2 P)(D 1 D 2 P) x = 

where a, b, ... , m are thirteen arbitrary bosonic parameters. In this paper, we focus on 
the simplest case where all of the parameters vanish: 

Ax =P t - P X (D X D 2 W) - (D X D 2 P)W X = 

(11) 

A 2 =W t - (D X D 2 W)W X + (-l)^ +1 A( J D lJ D 2 P)^ 2 ( J D 1 P)( J D 2 P)( J D lJ D 2 P) :[ . = 

We will subsequently refer to system (fTTT) as the supersymmetric polytropic gas dynamics 
equations. 

A real Grassmann algebra & is generated by a finite or infinite number of generators 
(£i)£2i • • •)• 1^ i s a graded vector space <& = <5 + <S X , where each variable Y in (3 that 
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we use is either even or odd. In general, for a Grassmann variable Y, we can define a 
"parity" Y which is if Y is even and 1 if Y is odd. In order to derive equation (11 II) . we 
have used the following Leibniz rule: 

de i (f9) = (d e J)g + (-iy(d ei g). (12) 

The partial derivatives with respect to odd coordinate satisfy the usual operational rules, 
namely doJOj = <5*. The operators d$., Qi, Q 2 , D\ and D 2 change the "parity" of the 
function acted on in the sense that it converts a bosonic function to a fermionic function 
and vice- versa. For example, d$.W is an odd superfield while de 1 dg 2 W is an even superfield, 
and so on. For further details, see the book by CornwelP^and the reference by DeWittP^ 



III Symmetries of the supersymmetric polytropic gas dynamics 
equation 



In order to determine the Lie superalgebra of infinitesimal symmetries of the supersym- 
metric equations (ITT]) , we make use of the apparatus of vector field prolongations as 
described in the book by P. J. Olver.^ Specifically, we make use of a generalization 
of this method to Grassmann variables in analogy with the case of the supersymmetric 
sine-Gordon equationPS 

A symmetry supergroup G of the system ffTTj) is a (local) supergroup of transforma- 
tions acting on the cartesian product of supermanifolds X x U, where X represents the 
independent variables (x,t,9i,9 2 ) an d U the dependent variables (superfields) (W,P). 
The action of G on the functions W(x, t, 9i, 9 2 ) and P(x, t, 9i, 9 2 ) maps solutions of ffTTT) 
to solutions of (jlip . Assuming that G is a Lie supergroup as described irpSHU] one can 
associate to it its Lie superalgebra of even left-invariant vector fields q, whose elements 
are the infinitesimal symmetries of the system ffTTT) . In particular, a local one-parameter 
subgroup of G consists of a family of transformations 

g £ : x l = X i (x,u,e), u a = U a (x,u,e) (13) 

where x = (x 1 , x 2 , x 3 , x 4 ) = (x, t, 9%, 62) are the independent variables and u = {u l ,u 2 ) = 
(W, P) the dependent variables. Here, e is a real parameter, whose range may be restricted 
depending on the values of x, t, 9 1 ,9 2 , W, P. Such a local subgroup is generated by a vector 
field of the form 

V = f(ll 4 + $(l ' M) i' (14) 
where the coefficients are given by 

C(x,u) = ^X% =0 , ^(x,u) = -^-U% =0 . (15) 
ae as 

The Lie superalgebra q, consisting of all vector fields of the form ( Fl4|) . then generates 
the Lie supergroup G. The advantage of working with the Lie superalgebra g instead 
of directly with the Lie supergroup G is that the equations defining the infinitesimal 
symmetries are linear. 

In order to determine the infinitesimal symmetries of a system of partial differential 
equations, it is useful to make use of the concept of the prolongation of a (super)group 
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action. The idea is that the transformation of the coordinates x l — > x\ u a — > u a induces 
a transformation of the derivatives 

-» <«» 

ox* ox* 

In order to make use of this concept, we define the multi-index J = (ji, ... , j p ), where 
= 0, 1, . . . and | J | = j i + . . . + j p . The space of coordinates on X x [7 is extended to the 
jet bundle Jk = {(x\u a ,Uj)\ |J| < £;}, which includes the coordinates and all derivatives 
of the dependent variables of order less than or equal to k. In our setting the jet bundle is 
a supermanifold as well, since X x U was. On the jet bundle, we define total derivatives 

A=^ + y"<7^-, (17) 
1 dx l ^ Jl du a j' y ' 

a, J J 

where J { = (j u . . . ,ji-xJi + hji+x, ■ ■ -,3n)- M ore generally, for J = (j u j 2 , ■ ■ ■ ,j n ), we 
define the composite derivative 

2? J = 2? 1 D 1 ---V, ■■■ V n V n ---V n . (18) 




The prolongation of a Lie (super)group action to the jet bundle Jk in turn induces a 
prolongation of the generating infinitesimal vector field in the Lie (super) algebra. For the 
vector field v given by (|14p . the k th order prolongation of the vector field v is 

pr«(v)=v+ £ <^(*y fc) )^, (19) 

where (f)j(x,u^ k ') are given by the formula 

^ = V j{r-C-Q-)+Cu% (20) 
or, equivalently, by the recursive formula 

The symmetry criterion (Theorem 2.31 in Olver^l) assumes that G is a connected Lie 
group of transformations acting locally onXx(J through the transformations 

Xi = X l (x,u,g), u a = U a (x,u,g), 

where g G G and A I/ (x, u^) is a non-degenerate system of partial differential equa- 
tions (meaning that the system is locally solvable and is of maximal rank at every point 
(xq, u^) G X x [/(")). Then G is a symmetry group of A = if and only if 

[pr (fc) (v)j (A) = whenever A = 0, (22) 

for each infinitesimal generator v of G. Using the results of one finds that the same 
criterion can be used also in the case of the Lie supergroup G and its Lie superalgebra of 
even left-invariant vector fields. 
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For the purpose of determining the Lie superalgebra of symmetries of the system ( TTTi) . 
let us write a vector field of the form 

v =£(ar, t, 6 1 ,6 2 , W, P)d x + r(x, t, 6 U 6 2 , W, P)d t + Pl (x, t, ls 6 2 , W, P)d Bl 

+ P2(x, t, 1; 2 , W, P)d d2 + A{x, t, 9 1 ,9 2 , W, P)d w + Q(x, t, 6 U 9 2 , W, P)d P , ' ' 

where £, r, A and Q are bosonic functions, while p\ and p 2 are fermionic. In general, we 
use the convention that the fermionic coefficients in a vector field expansion (in this case 
pi and p 2 ) precede the fermionic derivatives which they are multiplied by ( in this case 
dg 1 and de 2 respectively). 

According to the symmetry criterion, the vector field (1231) is an infinitesimal generator 
of the symmetry group of the system of differential equations fTTTj) if and only if 

pr^(v) [A k (x,t,e 1 ,9 2 ,W,P,W i2 \P^)] =0, k = l,2 (24) 

whenever A l (x,t,6 1 ,9 2 ,W, P,W^ , P^) = 0, I = 1,2. We use the following expressions 
for the total derivatives V x , V tl V 8l and V e , 2 : 

V x =d x + W x d w + P x d P + W xx d Wx + P xx d Px + W xt d Wt + P xt d Pt + W x9l d Wgi 
+ PxBxdp 9l + W X 2 dw e2 + P x e 2 d Pf)2 + W xxx d Wxx + P xxx d Pxx + W xxt d Wxt 
+ Pxxtd Pxt + W xx6l d Wx0i + P xxdl d PxBi + W xx e 2 d Wx02 + P xx e 2 d Pxf>2 + W xtt d Wtt (25) 
+ Pxudpu + W xtgi d Wtgi + P x te 1 d Pt g i + W xt e 2 d Wtf)2 + P x te 2 d Ptf)2 + W xdl d 2 d Weif)2 
+ Pxe^d^^, 

v t =d t + w t d w + p t d P + w xt d Wx + p xt d Px + w tt d Wt + p tt d Pt + w t9l d Wei 

+ Pte 1 dp Sl + W t e 2 dw B2 + P t e 2 d P02 + W xxt d Wxx + P xxt d Pxx + W xtt d Wxt 

+ Pxttdprt + W xWl dw xf)l + Pxteidp mei + W xt e 2 d Wx02 + P x te 2 d Pxf)2 + W m dwu (26) 

+ Ptttdp tt + W U e 1 dw t g 1 + -Pit0i<9p t9i + W t te 2 dw t e 2 + Ptte 2 dp t82 + Wto^dwe^ 

+ Pt6i6 2 dp gi92 , 

V Ql =d 6l + W 6l d w + P 01 d P + W x01 d Wx + P x6l d Px + W Wl d Wt + P tei d Pt + W e201 d Wg2 

+ Pe 2 Qxdpe 2 + W xx9l d Wxx + P xx e x d Pxx + W xtei d Wxt + P x te x d Pxt + W xd2dl d Wxg2 (27) 
+ Px0 2 eidp x e 2 + W tt o 1 dwtt + Ptte 1 dp tt + W t e 2 1 dw t93 + Pte 2 e-fip tH , 

and 

V e2 =de 2 + W d2 d w + P d2 d P + W xd2 d Wx + P x e 2 d Px + W t 9 2 dw t + P t e 2 d Pt + We 1 o 2 dw 6l 

+ P6x6 2 dp ei + W xx e 2 d Wxx + Pxx8 2 dp xx + W xt e 2 d Wxt + P x te 2 dp xt + W x9l e 2 dw xBl (28) 
+ Px6ie 2 dp x9i + Wue 2 dw t t + Ptte 2 dp tt + Wte^^Wt^ + Pt0i9 2 d Pt6i , 

We note that the chain rule for a Grassmann-valued composite function f(g(x)) i^2i|29| 

?L = ?9.?l. (29) 

The interchangeability of mixed derivatives (with proper respect to the ordering of odd 
variables) is of course assumed throughout. The second prolongation of the vector field 
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231) is given by 

pr (2) v =£ dx + rdt + + p2df)2 + + fia P + A 3 ^ + A'«9 m + A e ^d Wei 
+ A*^ + A xx d Wxx + A xt d Wxt + A^d Wx6i + A x92 d Wxe2 + A u d Wtt 

+ A^d Wt0i + A t9 >d Wte2 + A e ^d Weie2 + sra Pm + rtd Pt + n^d Pei (so) 
+ ^d PB2 + sr°d Pmm + n xt d Pxt + n^d Pxei + n^d PxS2 + n tt d Ptt 



+ n^d Ptei +n^d Pte2 +n e ^d P0i02 . 



Here, we use upper indices in coefficients A x , VL xdl etc. in order to distinguish them 
from partial derivatives, e.g. = de^Q. Applying the second prolongation (130]) to 

the equation ( II ip . we obtain the following conditions for the coefficients of the second 
prolongation f[3"Uj) 

Pl d Bl [P t - P X (D X D 2 W) - [D X D 2 P)W X ] 
+ p 2 d d2 [P t - P X {D X D 2 W) - [D X D 2 P)W X ] 
+ Q x d Px [P t - P X (D X D 2 W) - {D 1 D 2 P)W X ) 
+ tfd Pt [P t - P X (D X D 2 W) - [D X D 2 P)W X ] 
+ A xt d Wxt [P t - P X {D X D 2 W) - (D X D 2 P)W X ] 
+ A^dw t9i [P t - P X (D X D 2 W) - (D X D 2 P)W X ] 

+ A xd2 d Wx g 2 [Pt - P X (D X D 2 W) - (D X D 2 P)W X ] ^ 

+ A e ^d Wei02 [P t - P X (D X D 2 W) - (D X D 2 P)W X } 

+ n xt d Pxt [P t - P X (D X D 2 W) - [D X D 2 P)W X ] 

+ n t01 d Pt9i [P t - P X (D X D 2 W) - (D X D 2 P)W X ] 

+ n xe2 d Px02 [P t - P X (D X D 2 W) - [D X D 2 P)W X ] 

+ n ei92 d Peie2 [P t - P X (D X D 2 W) - (D X D 2 P)W X ] = 0, 



and 



Pl d 9l [W t - (D X D 2 W)W X + {-l)^ +1 A(D 1 D 2 Py- 2 (D 1 P)(D 2 P){D 1 D 2 P) x ] 
+ P2d d2 [W t - (D X D 2 W)W X + (-1)^ +1 A(D 1 D 2 P)^ 2 (D 1 P)(D 2 P)(D 1 D 2 P) X ] 
+ n x d Px [W t - {D X D 2 W)W X + (-iy +1 A(D x D 2 Py- 2 (D x P)(D 2 P)(D x D 2 P) x ] 
+ tfd Pt [W t - (D X D 2 W)W X + (-l)^ +1 A(D 1 D 2 F)T- 2 (D 1 P)( J D 2 P)(D lJ D 2 F) :l .] 
+ A xt d Wxt [W t - {D X D 2 W)W X + {-iy +1 A{D x D 2 Py<- 2 {D x P){D 2 P){D x D 2 P) x ] 
+ A t9 ^d Wtei [W t - (D X D 2 W)W X + (-1)^ +1 A(D 1 D 2 P)^ 2 (D 1 F)(D 2 P)(D 1 D 2 F) :E ] 
+ A xd2 d Wxe2 [W t - {D X D 2 W)W X + (-iy +1 A(D x D 2 Py- 2 (D x P)(D 2 P)(D x D 2 P) x ] 
+ A e ^d Wei62 [W t - (D X D 2 W)W X + (-l)T +1 A(D lJ D 2 P)^ 2 (D 1 P)(D 2 P)(D 1 D 2 P) a; ] 
+ Q xt d Pxt [W t - (D X D 2 W)W X + (-l)^ 1 A(D lJ D 2 P)^ 2 (D 1 F)( J D 2 P)(D 1 D 2 F) :E ] 
+ n^d PtBi [W t - (D X D 2 W)W X + (-iy +1 A(D x D 2 Py- 2 (D x P)(D 2 P)(D x D 2 P) x ] 
+ n xe *d Pxe2 [W t - (D X D 2 W)W X + {~iy +l A{D 1 D 2 Py- 2 {D 1 P){D 2 P){D 1 D 2 P) x \ 
+ n ^d P6i02 [W t - (D X D 2 W)W X + (-iy< +1 A(D x D 2 Py- 2 (D x P)(D 2 P)(D x D 2 P) x ] 
= 0. 

(32) 
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Note that proper respect to the ordering of fermionic terms is essential, e.g. A tdl is 
odd. We see that we only need to calculate the coefficients £l x , A xt , A tdl , A x02 , A 6 * 102 , 
Qxt ^ Qte 1 ^ Qxe 2 ^ QS1S2 j n equation ( 130]) . They are found from the superspace version of the 
formulae for the 1st and 2nd prolongation (|2Q|) 

A A = V A A -Y^V A ( B $ B , A AB = V B A A - Y,VbC C <$>ac, (33) 

B C 

where 

A,B,Ce {x, t, 9 Xi 6 2 }, C A = (£, r, p, a). (34) 

The derivation of these formulae is performed in the same way as in the bosonic case, 
working formally with infinitesimal transformations and keeping track of ordering prop- 
erties. The signs in the expressions below vary according to the considered equation 
Explicitly, the coefficients of the prolongation (1301) are given as follows: 

n x =a x + n w w x + n P p x - $ x p x - i w w x p x - £ p (p,) 2 - r x p t - r w w x p t 

- T P P x P t - pi, x Pe 1 - pi,wW x P$ 1 - pi,pP x Pe 1 - P2, x Pe 2 ~ P2,wW x P 02 (35) 

— P2,pPxPe 2 ; 



n* =n t + n w w t + n P p t - £ t P x - i w w t p x - £ P P x P t - T t P t - r w W t P t 

- T P {P t f - p 1>t P 9l - p 1>w W t P 6l - pi.pPtPe, - P2,tPe 2 - P2,wW t P 92 (36) 

- P2,pPtPe 2 , 



A xt =A xt + A xW W t + A xP P t + A tw W x + A ww W x W t + A WP W x P t + A w W xt 

+ A tP P x + A WP P x W t + A PP P x P t + A P P xt - UW X - i xW W x W t - i xP W x P t 

- txW xt - t m {W x y - Zww(W x ) 2 W t - iwp{W x ) 2 P t - 2£ w W x W xt 

- ttpw x p x - iwpw x w t p x - i PP w x p x p t - i P p x w xt - i P w x p xt - i t w xx 

- CwW t W xx - i P P t W xx - r xt W t - T xW (W t f - T xP W t P t - T x W tt - r tw W x W t 

- r ww (W t ) 2 W x - T WP W x W t Pt - 2T W W t W xt - r w W x W a - T tP W t P x 

- T WP (W t ) 2 P x - T P pW t P x P t - TpP x W tt - T P W t P xt - T t W xt - T P P t W xt 

- pi,xtW 8l - pi >xW W t W ei - Pi, x pP t W 9l - pi,twW x W 6l - pi, x W tei - pi,tW xdl 

- pi,wwW x W t W dl - pi, W pW x W 6l P t - pi,wW xt W 9l - pi, w W Wl W x 

- Pi,wW x6l W t - pi, tP P x W 6l - pi, W pPxW t W ei - p lt p P P x P t W ei - pi t pP xt W ei 

- Pl,pPxW t9l - p 1 , P PtW x g 1 - P2,xtWg 2 - p 2 ,xwW t Wg 2 - p 2 , x pPtWg 2 

~ p2,twW x Wg 2 - P2,xW t g 2 - p 2 ,tW x 2 - P2,WwW x W t Wg 2 - P2,WpPtW x Wg 2 
~ p2,wW xt Wg 2 - p2, W W W2 W x - p2, W W x e 2 W t - p 2 ,tpPxWg 2 - p2, W pPxW t Wg 2 
~ P2,PpPxPtWg 2 - P2.pPxtWg 2 - p2,pW t g 2 P x ~ p2,pW x g 2 P t , 

(37) 
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=A t01 + A tw W ei + A tP P 01 + Ag lW W t + A ww W t W ei + A WP W t P ei + AwW t(?1 
+ A eiP P t + A WP P t W 6l + A PP P t P 9l + A P P t9l - i m W x - CtwW x W 6l 

- ZtpW x P 01 - £ t W x01 - Ce lW W x W t - t ww W x W t W 6l - iwpW x W t P 9l 

- ZwW t W x9l - ZwW x W tei - £ ei pW x P t - twpPtW x W 9l - tppW x P t P 9l 

- tpPtW x9x - tpW x P t6x - ie x W xt - iwW xt W 9l - &W xt P 01 - r Wl W t 

- r tw W t W 9l - r tP W t P 9l - r t W t9l - r 9lW (W t ) 2 - r ww (W t ) 2 W 9l 

- r WP {W t ) 2 P ei - 2r w W t W m - r 9lP W t P t - r WP P t W t W dl - T PP W t P t P 9l 

- T P P t W t9l - T P W t P t9l - T 9l W tt - T W W tt W 9l - T P W tt P 9l - p ht01 W ei 
+ Pi,tpPe 1 W 9l - pi, 9lW W t W 9l + pi,wpPe 1 W t W ei + pi,wW t01 W 01 

- pi,e 1 pPtW 9l + pi,ppP t Pe 1 We 1 + pi,pP t e 1 W 01 - pi,e 1 W t01 + pi,wW 01 W t01 
+ pi,pPe 1 W Wl - p 2 ,te 1 Wo 2 + p2,twW 9l W 92 + p 2 ,tpPe 1 W 92 - p 2 ,tW 0102 

- p2,e 1 wW t W 92 + p2,wwW t W 9l W 92 + p 2 ,wpPe 1 W t Wg 2 + p 2 ,wW t01 W 02 

- p2,wW t W 9l92 - p2, ei pPtW 92 + p2,wpPtW 9l W 02 + p2,ppP t Pe 1 W 92 

+ p2,pPte 1 W 92 - p2,pPtW 9l92 - p2,e 1 W t92 + p2,w w di W te 2 + p2,pPe 1 W t92 , 



A x92 + A xW W 92 + A xP P 92 + A 92W W X + A WW W X W 92 + A WP W X P 02 + A W W X02 
+ A 92P P X + A WP P X W 92 + A PP P X P 92 + A P P x92 - £ X02 W X - i xW W x W 92 - £ xP W x 

- txW x02 - i 92W {W x f - iww{W x ) 2 W 92 - iwp{W x ) 2 P 92 - 2£ W W X W X02 

- ie 2 pW x P x - iwpP x W x W 92 - i PP W x P x P 92 - t P P x W x92 - i P W x P x92 

- ie 2 W xx - £ W W XX W 02 - i P W xx P 02 - r x92 W t - r xW W t W 92 - r xP W t P 92 

- t x W W2 - T 9 , 2W W x W t - T WW W x W t W 92 - T WP W x W t P 92 - r w W x W t92 

- r w W t W x92 - r d2P W t P x - T WP P x W t W 92 - r PP W t P x P 02 - r P P x W t92 

- r P W t P x02 - T 02 W xt - T W W xt W 02 - T P W xt P 02 - pi, x02 W 0l + pi, xW W 02 W 01 
+ pi, x pPe 2 W 01 + pi, x W 0102 - pi, 02 wW x W 01 + pi,wwW x W 02 W 01 

+ pi,wpPe 2 W x W 01 + pi, w W x02 W 01 + pi, w W x W 0102 - pi, 02 pP x W 01 

+ Pi,wpPxW 02 W 01 + pi, PP P x P 02 W 01 + pi, P P x02 W 01 + pi, P P x W 0l02 - pi, 02 W x01 

+ Pi,wW 02 W x6l + Pi, P P 02 W x01 - p2, x e 2 W 02 + p 2 ,xpPe 2 W 02 - P2,e 2 wW x W 02 

+ p2,wpPe 2 W x W 02 + p2,wW x02 W 02 - p2,e 2 pP x W 02 + p2,ppP x Pe 2 W 02 

+ p2,pP x e 2 W 02 - p 2 ,e 2 W x02 + p2,wW 02 W x02 + p 2 ,pPe 2 W x02 , 
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=A 0102 - A 0lW W 02 - A ei pPe 2 + A 02W W 01 - A WW W 01 W 02 - A WP W 01 P 02 
+ A W W 0102 + A 02 pP 01 — A WP P dl Wg 2 — A PP P 01 P 02 + A P P 0102 — £ 0102 W X 
+ ie lW W x W d2 + ie lP W x P d2 + i 01 W x02 - Ce 2 wW x W 01 + ZwwW x W 01 W 02 
+ CwpW x W 6l Pe 2 + ZwW ei W x e 2 - ZwW x W ei0a - & 2 pW x P 6l 
+ iwpPe x W x W d2 + ip P W x P ei Pe 2 + ZpP dl W x02 - ^ P W X P 0102 - £e 2 W x01 

- iwW d2 W xdl - i P Pe 2 W xdl - r 0102 W t + r 0lW W t W 02 + r 0lP W t P 02 + T ei W m 

- T0 2W W t W ei + r ww W t W 01 W 02 + T WP W t W 01 P 02 + T W W 01 W t02 

- T W W t W eid2 - T0 2 pWtP ei + T W pP ei W t We 2 + rppW t P ei Pe 2 + r P P 01 W t02 

- T P W t P ei e 2 - Te 2 W tdl - r w W 02 W t01 - T P P 02 W t01 - pi, 0102 W 01 
+ Pi,e 1 wW 9l We 2 + pi, 01 pW 0l P 02 - p\^W ^ 2 + p\ )02 pP 0x W 01 

- pi,wpPe 1 W ei W d2 + pi,PpPe 1 Pe 2 W 6l - pi, P W 01 P 0102 - p hP P 01 W 0102 

- p2,e 1 e 2 Wg 2 - p 2 ,e 1 pPe 2 Wg 2 + p2,e 2 wW 01 W 02 - p2,wpPe 2 W 01 W 02 

+ p2,e 2 pPe 1 We 2 + p 2 ,ppPe 1 Pe 2 W 02 - p 2 ,pW 02 P 0102 - p 2 ,e 2 W 0102 + p2,pPe 2 W 0102 , 

(40) 



&xt + &<xpPt + QxwWt + Vt t pP x + VtppP x Pt + Vtpy/PxWt + VtpP xt 

+ n tw w x + n PW w x p t + n ww w x w t + n w w xt - i xt p x - i xP p x p t - ixwP x w t 

- t x Pxt - itp(Px) 2 - ip P {Px?Pt - ip W {Px?w t - ii P p x p xt 

- itwPxW x - i PW PxPtW x - twwPxW x w t - i w w x p xt - i w p x w xt - i t p xx 

- tpPtPxx - twW t P xx - r xt P t - T xP {P t f - r xW P t W t - T x P tt - r tP P x P t 

- T PP (Pt) 2 P x - r PW P x P t W t - 2r P P t P xt - r P P x P tt - r tw P t W x 

- T PW {P t ) 2 W x - r ww P t W x W t - r w W x P tt - r w P t W xt - r t P xt - T W W t P xt 

- pi,xtPe! — P\,xpPtPe x — Pi,xwW t P 01 — PijpPxPe! — P\,xPte x — p\,tPxe x 

- p\,ppPxPtPex — Pi,pwPxPeiW t — p\,pP x tPe x — p\,pPte x Px 

- p\,pPxd x Pt - Pi,twW x P 01 - pi,pwW x P t P 01 - pi,wwW x W t P 01 - pi,wW xt P 01 

- Pi,wW x P t01 — pi : \yW t P x01 — P2,xtPe 2 — p2,xpPtPe 2 — p2,x\yW t P 02 

- p2,tpP x Pe 2 — P2, x Pte 2 ~ P2,tPx8 2 — p2,ppP x PtPe 2 — p2,pwW t P x Pe 2 

- p2,pPxtPe 2 — P2,pPte 2 Px — P2,pPxe 2 Pt — P2,twW x P 02 — P2,pwW x P t P 02 

- p2,WwW x W t P 02 - p2,wW xt P 02 - p2,wPtd 2 W x - P2,wPx9 2 W t , 

(41) 
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=n t6l + n t pPe 1 + to t wW ei + n eiP p t + n PP p t p ei + n PW p t w ei + n P p tei 
+ n 6lW w t + n PW w t p 6l + n w w t w ei + n w w Wl - i m p x - i tP p x Pg, 

- itwP x We x — itPxdx — ie x pPxPt — ippPxPtPe x — ipwPxPiWe x 

- ipPtPxe, - ipPxPte, - fawPxWt - ti PW W t P x P dl - iwwPxW t W 8l 

- ZwW t P x6l - twPxW tQl - i H P xt - i P P xt P dl - iwPxtWe, - r m P t 

- r tP P t P dl - T tw P t W ei - T t P m - r eiP (P t ) 2 - T PP {P t fP 6l 

- r PW (P t ) 2 W ei - 2r P P t P tdl - T 9lW P t W t - r PW W t P t P 6l - r ww P t W t W dl 

- T W W t P t6l - T W P t W tei - T 6l P tt - T P P tt P 6l - T W P tt W dl - Pi^Po! 

+ Pi,twWg 1 Pg 1 — p\fi xP PtPe x + Pi,pwW ei P t P 9l + pi t pPte 1 Pe 1 

- Pifi 1 wW t Pe 1 + Pi,wwW t Wg 1 Pg 1 + Pi,wW t9l Pg 1 - pi,g 1 P t e 1 + Pi,pPe 1 Pte 1 
+ Pi,w^6\Pte^ — P2,t6iPe 2 + p2,tpPe 1 Pe 2 + p2,twW 01 Pg 2 — p2,tPe 1 8 2 

- P2,e 1 pPtPe 2 + p2,ppPtPe 1 Pe 2 + p2,p\vWg 1 P t Pg 2 + P2,pPte 1 Pe 2 

- P2,pPtPe 1 e 2 ~~ P2,o 1 wW t Pe 2 + P2,p\vW t Pg 1 Pg 2 + p2,wwW t Wg 1 Pg 2 

+ p2,wW t e 1 Pe 2 — P2,wW t Po 1 o 2 — p2,e 1 Pte 2 + P2,pPe 1 Pte 2 + p2,wWg 1 P t g 2 , 



—Qx9 2 + &xpPe 2 + Q x wWo 2 + Qg 2P P x + Vt PP P x Pg 2 + VL P wP x Wg 2 + Q P P x g 2 

+ ^e 2 wW x + QpwW x Pg 2 + Vt ww W x Wg 2 + Vl w W x g 2 — £ x e 2 P x — i xP P x Pe 2 — £, x wPxWg 2 

- ixP x e 2 — ie 2 p{Px) 2 — ipp{Px) 2 Pe 2 — £,Pw(Px) 2 Wg 2 — 2£ P P x P x g 2 

- ie 2 wPxW x - £ P wW x P x Pe 2 - iwwPxW x We 2 - £ w W x P x e 2 - i w PxW x g 2 

- ie 2 Pxx — ipPx X Pe 2 — £wPxxW e2 — r x g 2 P t — T xP P t P e , 2 — T xW P t W d2 

- T x P t g 2 - Tg 2 pP x P t - T PP P x P t Pe 2 ~ T PW P x P t Wg 2 - T P P x P t e 2 

- T P P t P x02 ~ Tg 2W P t W x - T PW W x P t Pg 2 - T WW P t W x Wg 2 - T W W x P t e 2 

- T W P t W x e 2 - Te 2 P xt - r P P xt Pg 2 - T W P xt We 2 - pi, x e 2 Pe 1 + Pi, x pPe 2 Pe 1 
+ Pi,x\vWe 2 Pe 1 + P\, x Pe^ 2 - pi,e 2P P x Pe 1 + Pi,ppP x Pe 2 Pe 1 

+ Pi,pwWg 2 P x Pg 1 + P\^ p P x q 2 Pq x + Pi,pP x Pe!0 2 — Pi,e 2 wW x P 6l 

+ pi,pwW x Pg 2 P 8l + pi,wwW x Wg 2 P 6l + pi,wW x g 2 P ei + pi,wW x P 6l e 2 - pi,e 2 P x e 1 

+ Pl,pP0 2 PxOi + Pl,wWg 2 P x g l — P2,x9 2 Pe 2 + p2,xwWg 2 Pg 2 — p2fi 2 pP x Pe 2 
+ P2,PwWg 2 P x Pg 2 + p2,pPx6 2 Pe 2 ~ p2,e 2 wW x Pg 2 + p2,WwW x Wg 2 Pg 2 
+ P2,wW x g 2 Pg 2 — p2,e 2 Px6 2 + P2,pPe 2 Px9 2 + p2,wWg 2 P x g 2 , 

(43) 
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£l dl02 =£Iq i62 — Vtg lP Pg 2 — £lg lW Wg 2 + Q,0 2 pPe 1 — £lppPe 1 Pe 2 — ^pwPb^b 2 

+ QpPe x B 2 + — ^PwWg 1 P$ 2 — Vt ww Wg 1 Wg 2 + QwWg^ — ig^Px 

+ ie±pPxPo 2 + £eiwPxWo 2 + £g 1 P x e 2 ~ ^e 2 pPxPe x + £,ppPxPe 1 Pe 2 

+ ipv/PxPe 1 We 2 + ipPe 1 P x e 2 — t,pPxPe x e 2 — io 2 wPxW Bl 

+ £,PwWo 1 P x Pe 2 + iwwPxWe 1 We 2 + ^,wWg 1 P x g 2 — ^wPxWg ± g 2 — £g 2 P x e 1 

~ £pPe 2 Px8i — £,wWg 2 P x g 1 — Tg x g 2 Pt + Tg lP P t Pg 2 + Tg lW P t Wg 2 + T 8l P t g 2 

— Te 2 pP t P ex + T P pP t P Bl Pe 2 + rjwiVWfc + TpP dl P t g 2 

— rpP t P ei g 2 - T0 2W P t W dl + Tp W W ei P t Pe 2 + r w P t W ei W e2 + r^W^P^ 

— T wPt^e 1 e 2 — T d 2 Ptd x — T pPe 2 Pte 1 — T wWg 2 P t g 1 — p± ) g 1 g 2 Pg 1 
+ Pi,e 1 pPe 1 Pe 2 + Pi^w-Pi^ 7 ^ — pi^Pe^ + Pi,e 2 wWg 1 Pg 1 

— Pl,P\vWg 1 Pg 1 Pg 2 + Pl,WwWg 1 Wg 2 Pg 1 — px t w Pqi Wg x g 2 — Pi ) wWg 1 Pg 1 g 2 

— P2,e 1 e 2 Pe 2 — P2,9 1 wWg 2 Pg 2 + p 2> g 2 pPg 1 Pg 2 — P2,P\vWg 2 Pg 1 Pg 2 

+ P2.8 2 wWg 1 Pg 2 + P2,WwWg 1 Wg 2 Pg 2 — P2,wPe 2 Wg 1 g 2 — p 2 fi 2 Pg x g 2 + P2,wWg 2 Pg 1 g 2 , 

(44) 



Substituting the above formulae into the equations (T5T|) and (1521) and making the 
substitutions from the original supersymmetric equations ffTTl) 

P t = P X {D X D 2 W) + (ADaP)^, 

(45) 

W t = (D 1 D 2 W)W X + (-iyA(D 1 D 2 Py- 2 (D 1 P)(D 2 P)(D 1 D 2 P) x , 

we obtain a series of determining equations for the functions £, t, pi, p 2 , A and in (1231) . 
We provide a specific solution of these determining equations, with the form: 

X ) = 2<7n + C 2 - r{t, 9 2 ) = 2C 3 t + C 4 - Krft, 

pM = + Ki, p 2 (6 2 ) = C 3 9 2 + K*, A(W) = (3d - C 3 )W, 

where Ci, C 2 , C 3 are bosonic constants, while K\ and K 2 are fermionic constants. Here, 
the function Q(P, 7) is parametrized by the polytropic exponent 7. Thus, the superalgebra 
£ is spanned by the following vector fields: 



P 1 = d x , P 2 = d t , Qi = -0 1 d x + d 0l , 
L x (7) = 2xd x + 9 1 dg 1 + 3Wd w + (gf ) Pd P , 

l 2 ( 7 ) = 2^ + e 2 d d2 - wd w + (3=3) pd P . 

Here, the generators Pi and P 2 represent translations in the x and t directions respectively, 
£1(7) and 1/2(7) are two families of dilations in the independent and dependent variables, 
and Qi and Q 2 are the supersymmetric transformations already determined in equation 
(EJ. The dependence on 7 is reflected in the form of the two dilations, £1(7) and L 2 (-y), 



Q 2 = -6 2 d t + dg 2 , 

(47) 
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Table I: Supercommutation table for the Lie superalgebra £ spanned by the vector fields (|4"T|) 





Li 


Pi 


Qi 


L 2 


p 2 


Q 2 


Li 





-2Pi 


-Qi 











Pi 


2Pi 

















Qi 


Qi 





-2Pi 











L 2 














-2P 2 


-Qi 


P 2 











2P 2 








Q 2 
















-2P 2 



which we abbreviate for convenience as L\ and L 2 respectively. The supercommutation 
relations of the superalgebra generators (j47p are given in Table I. 

This superalgebra differs in structure from a classical Lie algebra in the sense that the 
diagonal contains non-zero elements, i.e. {Qi,Qi} = —2Pi and {(3 2 ,Q 2 } = — 2P 2 . 

We now proceed to classify the one-dimensional subalgebras of the Lie superalgebra 
into conjugacy classes according to the action by the Lie supergroup generated by the 
Lie superalgebra £. We construct a list of representative subalgebras such that each 
subalgebra of £ is conjugate to one and only one element of the list. We perform this 
analysis using the methods described in (for example).®^ We begin our classification by 
writing the superalgebra £ as a direct sum of two semi-direct sums of smaller algebras in 
the following way 

£ = {{LJ 35 {P 1 ,Q 1 }} © {{L 2 } 35 {P 2 ,Q 2 }}, (48) 

Consider first the semi-direct sum 

a = {L 1 }^{P 1 ,Q 1 }. (49) 

This algebra a can be classified by using the method for semi-direct sums of algebras.^ 
This method leads to splitting subalgebras (of the form m D n, where m C {^i} and 
n C {Pi, Qi}) and non-splitting subalgebras which cannot be written in this form. 

Let us first discuss splitting subalgebras. The subalgebras of {£1} are {0} and 
For each of these two subalgebras, we must consider all subalgebras of {Pi, Qi} which are 
invariant under (super) commutation under the given subalgebra of {Li}. 

For subalgebra {0}, all subspaces of {Pi, Qi} are invariant subalgebras. That is, 

[aPi + /xQi,0] = C {0} (50) 
So the possible representative subalgebras are {0}, {-Pi}, {fiQi}, {Pi + ^Qi} and {Pi, Qi} 

We classify such subalgebras in the following way. We begin with the trivial subalgebra 
a o = {0}. Let us now consider the one- dimensional subalgebra {Pi} and apply the action 
of the subgroup e^ Lx,Px ^ 1 ^ generated by a to it. If we take any arbitrary element of a: 

Y = aLx + pPx + ijQi, (51) 
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then 

[Y,P 1 ] = -2aP 1 . (52) 

Therefore, under the action of the subgroup e^ Ll ' Pl ' Ql \ we obtain Pi — > e~ 2a P\. So, 
Oi = {-Pi} is the only subalgebra in its conjugation class. 

Consider now the subalgebra {nQ\}. If Y = aL\ + (3 Pi + T]Qi, then we obtain 

[Y,fxQx] = -a^Qi + rjjiPi. (53) 
Consequently, under the group action of e^ Ll ' Pl ' < ^ 1 ^ , we get 

if a * ° (54) 
V ; 



~ \ fiQi + wPii if a = 



So, by re-scaling the rj, we can obtain other multiples of the generator fiQi + i][iP\. We 
represent all these possibilities by the subalgebra a 2 = {^Qi}- 

The subalgebra {Pi + fiQi}, under the action of the group, leads to 

{Px + fxQi} -> {/xQi + (1 + w)Pi}- (55) 
All of these are covered under the subalgebra a 3 = {Pi + fiQi}. 

The remaining subalgebra, {Pi,<5i} is already two-dimensional and is therefore not 
included in our classification of one-dimensional subalgebras of a. 

The subalgebra {£1} is already one-dimensional, so the only splitting one-dimensional 
subalgebra of a obtained from consideration of the subalgebras of {Pi, Qi} invariant under 
{Pi} is a 4 = {Li} itself. 

Let us now discuss non-splitting subalgebras of a. We consider vector spaces of the 
form 

V = {L 1 + a 1 P 1 + a 1 Q 1 }. (56) 
Applying a cocycle to {Pi}, we obtain 

U + Ai[Pi, la] + Aa[Qi, £1] = Li + 2X 1 P 1 + A2Q1 (57) 

so, if we let A x = —\a\ and A2 = — a 2 , then for (J56"|) . we get 

V -> {Li + 2X X P X + AsQi + c x P x + AQi} = {Pi} (58) 

Thus, there are no non-splitting subalgebras of o. So, the algebra a = {Pi, Pi, Qi} has 
the one-dimensional subalgebra classification 

a = {0}, ai = {Pi}, o 3 = {/xQ 1 }, a 3 = {Pi + M<2i}, a 4 = {Pi} (59) 

Similarity, the algebra b = {P 2 ,P 2 ,<3 2 } has the one-dimensional subalgebra classification 

b = {0}, bi = {P 2 }, b 2 = {^Q 2 }, b 3 = {P 2 + //Q 2 }, b 4 = {P 2 }. (60) 

Considering now the direct sum 

£ = a e 6 = {Pi, Pi, Qx} e {p 2 , p 2 , Q 2 } (ei) 
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we use the Goursat method for the classification of direct sumsP^"^ 



The one- dimensional non-twisted subalgebras are formulated by successively taking 
the direct sums of each subalgebra of a with each subalgebra of b. This results in the 
following list: 

£i = {Pi}, £ 2 = {M2i}, £ 3 = {Pi + /iQi}, £ 4 = {Pi}, m 
£ 5 = {P 2 }, £ 6 = kQ 2 }, £7 = ^2 + ^2}, £s = {L 2 }, ((> " 



For twisted subalgebras, there are 16 possibilities of combining two subalgebras to- 
gether. Two subalgebras, Oj and bj can be twisted together if a homomorphism r exists 
from dj to bj\ 

r(Oi) = bj. (63) 

For example, if we consider Oi = {-Pi} and b% = {P 2 }, then the most general homomor- 
phism from di to bi is: 

t:P 1 ^ kP 2 , (64) 

where k is any scalar. This gives us the family of twisted subalgebras V = {Pi + kP 2 }. 
By considering the action of the full supergroup e^ Ll ' L2 ' Pl ' P2 ' ( ^ 1 ' ( ^ 2 ^ on V, we see that k 
can be rescaled by any positive real number. So, the representative subalgebra is: 

£ 9 = {P 1 + £ P 2 }, e = ±l. (65) 

Considering all possibilities, we obtain the list of all representative one-dimensional sub- 
algebras of the superalgebra £ generated by the vector fields (J4] 
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£1 = {A}, 


£2 = 


= {^Qi}, 






£ 3 = {Pi + ^Qi}, 


£ 4 = 


= {Li}, 






£5 = {P2}, 


£ 6 = 


= {KQ2}, 






£ 7 = {P 2 + z/Q 2 }, 


£ 8 = 


= {L 2 }, 




£9: 


= {P 1 + eP 2 ,e = ±l}, 


£10 


= {Pl+KQ2}, 


£11 


= {PI- 


^sP 2 + uQ 2 ,£ = ±1}, 


£12 


= {L 2 + eP 1 ,e = ±l}, 






£l3 = {P2 + //Q1}, 


£14 


= {[K>i +KQ2}, 




£15 


= {P 2 + {±Qi + vQ 2 }, 


£l6 


= {L 2 


£17 


= {P1- 


\-eP 2 + (iQi } e = ±1}, 


£l8 


= {Pi+ {1Q1 + KQ2}; 


{Pi 


+ eP 2 4 


- nQi + uQ 2 ,e = ±1}, 


£20 


= {L 2 + eP x + (jQi,e = ±1} 




£21 : 


= {L 1 +eP 2 ,e = ±l}, 


£22 


= {Li +KQ2}, 


£23 


= {Li- 


f eP 2 + vQ 2 ,e = ±1}, 


£24 


= {Lt + k^k^O}. 



This list of subalgebras will allow us to systematically use the symmetry reduction method.^ 



IV Symmetry reductions and invariant solutions 

We now make use of the symmetry reduction method in order to find invariant solutions 
of the supersymmetric polytropic gas dynamics equations ffTTl) . Replacing the specific 
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Table II: Invariants and change of variable for subalgebras of the Lie superalgebra £ spanned 
by the vector fields ([37]) 



Subalgebra 


Invariants 


Superfields 


£1 = mi 


t, 9 X , 9 2 , W, P 


W = W(t,9 1 ,9 2 ),P = P(t,9 1 ,9 2 ) 


£ 3 = {Pi+M2i} 


t, n = e 1 - nx, e 2 , w, p 


W = W(t,n,9 2 ),P = P(t,n,9 2 ) 


£ 4 = {Lt} 


t, t 1 = x-V 2 9 1 , 9 2 , 

( T + 5 A 

A = x~ 3/2 W, B = x \^+ 2 JP 


W = x*l 2 A{t,T 1 ,9 2 ), 
( 7+5 

P = aA 2 T+ 2 /£(i, 71,02) 


£5 = {^2} 


x, 9 U 9 2 , W, P 


W = W(x,0i,0 2 ), P = P{x,9 x ,9 2 ) 


2 7 = {P 2 + vQ 2 } 


x, 9 x ,r 2 = 9 2 -vt, W,P 


W = W(x,9 1 ,t 2 ),P = P(x,9 1 ,t 2 ) 


£ 8 = {L 2 } 


x, 9t,r 2 = t~ 1 l 2 9 2 , 


W = t- l l' 2 A{x,9 1 ,T 2 ), 




(1-3) 

A = t 1/2 W, B = x \ 2 ~<+ 2 )p 


( ~<-' A \ 

P = x\ 2 ~< +2 )B{t,9i,T 2 ) 




a = x-et, 6»i, 9 2 , W, P 


W = W(a,9 1 ,9 2 ),P = P(a,9 1 ,9 2 ) 


£10 = {Pi +EQ2} 


a = t + v9 2 x, Ui, 
T2 = 9 2 - v_x, W, P 


W = W{a,9i,T 2 ), P = P(a,9i,T 2 ) 


£11 = {Pi + eP 2 + v_Q 2 } 


a = t — ex + v_x9 2 , 9\, 

f\ TT7" 7~) 

t 2 = 9 2 - vx, W , P 


W = W{a,9x,r 2 ), P = P((7,0i,r 2 ) 


^12 ~~ lP2 t j 


n — T — iplnt #1 To — t~ 1 / 2 9n 


vv — i »Ai^y,(7i, i 2 j 

P = t(^)B(a, 9i,t 2 ) 


£13 = {P2 + MQ1} 


a = x + n = 9\ — /it, 


W = W(o,t 1 ,9 2 ) 

P — Pin- t, 


£ie = {L 2 + ±±Qi} 


= x + \fi9\ lni, ri = 9\ — \fi\wt, 


W = t- 1 ' 2 A{a,T 1 ,T 2 ) 




T 2 — I ' (72, J\ — t 1 VV , 


-t — C v y t) I^CT, Tl , T2 ) 










<T — FT — / -1- lltf)i 

n = 0i - e/ut, C W, P 


W — W((T T1 9<>) 

P = P((T,ri,0 2 ) 


£20 = {^2 +ePi +vQi} 


cr = x + ^/U^i hit — ielni, 


= r 1 / 2 < 4( < 7,r 1) T 2 ) 




t\ = 9\ — \ filnt, 
T2 = t~ l / 2 9 2 , A = t l / 2 W, 


P = t(^) J B(a,ri,r 2 ) 




( 7—3 ^ 

B = t~\*T+*)p 




£21 = {^1 + eP 2 } 


a = t — ^elnx, t\ = x~ x l 2 9\, 


M/ = x 3 / 2 ^(a,ri,0 2 ), 




e 2 , a = x-v 2 w, 


P = xV^^( f j,Ti,0 2 ) 




B = x~vw*)p 




£22 = ^1+^2} 


a = t+ \v_9 2 \nx, T\ = x~ 1/2 9\, 


W = x 3 /M(a,ri,r 2 ) 




t 2 = 9 2 - \u\nx, A = x~V 2 W, 


P = xV^+ 2 /0(cr,ri,T2) 




B = x"(^)p 




£23 = {il+£P 2 + l£Q2} 


a = t + \v9 2 In x — \e In x, 


VF = x 3 /M( ( J,r 1 ,r2) 




n = x -i/2# 1; 
t 2 = 9 2 - \u\nx, A = x~ 3 / 2 W, 

B = z"(^)p 


P — nr\ 27 + 2 } Ipf ifj 7--, t-,,"! 






£24 = {Li + /cL 2 } 


1 j_ 

a = t kx, T\ = t 2k 9\, 

1 k— 3 
T 2 = t"26» 2 , A = t~2kW, 

/( 7 +5)+fc( 7 -3)\ 
B = t \ 2fc( 7 +l) ) p 


'i—k 

W = t~2k~ A(a,n,T 2 ) 

/( 7 +5)+fc( 7 -3)\ 
P = A ^(a,Tl,T 2 ) 
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forms of the superfields W and P as listed in Table II into the supersymmetric equation 
pip we obtain symmetry reductions. 

In most cases these reductions involve heavy computations and we omit them for this 
reason. So, we present in this paper only certain classes of the invariant solutions involving 
(in most cases) freedom of arbitrary functions of one or two variables whose arguments 
can be either bosonic or fermionic. We propose two types of solutions. The first type has 
the form 

w = s(t F(e 1 , e 2 )) , p = R(t G(e u e 2 )) , (66) 

where F and G are arbitrary functions of the fermionic arguments 9± and 9 2 and £ is the 
symmetry variable. In some cases, we specify them to be in the form 

F(9 1 ,9 2 ) = aF 1 (9 1 ) + PF 2 (9 2 ) and G(0 U 6 2 ) = «G?i(0i) + \G 2 (9 2 ), (67) 

where a, (3, k and A are fermionic constants. The first type of solution is of the form 

W = S& F(9 U 9 2 )) , P = R (£, G(9 U 9 2 )) , (68) 

and corresponds to the results which we have obtained for subalgebras £i, £3, £7, £g, 

£l0) £ll, £l2 ; £l3) £l6) £l7) £20) £21; £23 an d £24- 

The second class of solution is of the type 

W = S^,9 U 9 2 ), P = R(£,9 1 ,9 2 ), (69) 

which we develop as a Taylor expansion with respect to 9\ and 9 2 with non-constant 
coefficients which depend on the symmetry variable £. Since all terms that include (#i) 2 
or (9 2 ) 2 vanish, the solutions (1691) are of the form 

W = Mi(0 + 9 1 M 2 (0 + 9 2 M 3 (0 + 0i0 2 M 4 (O, (7) 
P = N 1 (£) + e 1 N- l (Q + e 2 N 3 (Z) + 9 1 6 2 N i (Q [ j 

Replacing into the reduced equations, we obtain the solutions corresponding to subalge- 
bras £4, £ 5 , £ 9 , and £ 2 2- 

We obtain four different general classes of solutions. 



A Solutions involving only fermionic variables 

For the subalgebra £1 = {-Pi}, the reduced system of equations is 

P t = 0, W t = 0, (71) 

and the obtained solution is 

W = W(9 1 ,9 2 ), P = P(9 1 ,9 2 ). (72) 

This solution has two arbitrary functions of the two fermionic variables 9\ and 9 2 . Using 
the Taylor expansion of ( 1721 in 9\ and 9 2 , the solution ( 1721 can be expressed as 

W(p x , 9 2 ) = Ci + Crfx + C39 2 + CM, 

(73) 

P(9i, 9 2 ) = C 5 + C 6 0! + C 7 9 2 + C 8 9 1 9 2 , 



18 



where Ci, C4, C5, Cs are bosonic constants, and C2, C3, C*6, C7 are fermionic constants. 

For the subalgebra £5 = {P2}, the reduced system of equations is 
P x W 6l e 2 - 9 x P x W x e 2 + P 8l e 2 W x - d x P^W x = 0, 

W 0102 W X - e 1 W x e 2 W x + (-iyA(-P 9l e 2 y- 2 [P 01 Pg 2 P xei e 2 + P 9l Pe 2 0iP xx e 2 + 9 1 P x Pe 2 P x01 e 2 ) 
- Afr - 2)9 1 (P ei e 2 y- 3 P x e 2 Pe 1 Pe 2 P x e 1 9 2 = 

(74) 

and the obtained solution in the fermionic sector is 

W = gUm - fc)0 x + (// 2 - fa)B*), P = f(W) = /((//! - + (fi 2 - f3 2 )6 2 ), (75) 

This solution has two arbitrary functions, / and g, of one fermionic argument (/ii — 0i)9i + 
(/i 2 — 2 )9 2 . Using the Taylor expansion, the solution f!75|) can be reduced to the form 

W = C + Ci [(/ix - /3i)0i + - fo)6 2 ] , P = 02 + 03 [(//! - 000! + (0 2 - fr)6 2 ] , 

(76) 

where Co and C 2 are bosonic constants and C\ and C3 are fermionic constants. 

For the subalgebra £9 = {Pi + ei^}, the obtained solution, 

W = fi-axOx - 0262), P = gi-axOx - 0262), (77) 

depends only on the fermionic sector 9\ and 62- This solution has two arbitrary functions, 
/ and g, of one bosonic argument —ol\9\ — 026*2. In analogy with the case £5, this solution 
can be reduced to 

W = Co + Ci [-Q10! - 02^2] , P = C 2 + C A [-O10! - «20 2 ] , (78) 
where Co and C2 are bosonic constants and C\ and C3 are fermionic constants. 

For the subalgebra £17 = {Pi + 6P2 + l^Qi}, the obtained solution in the fermionic 
sector is 

W = f(9 2 )+(j.(6 x ), P = g(9 2 )+ [ x(9 1 ) (79) 

which contains two arbitrary functions, / and g, of one fermionic argument 9 2 . This 
solution can be written in the form 

W = C + Cj0 2 + //ft, P = C 2 + Cs9 2 + ^1, (80) 

where Co and C2 are bosonic constants and C\ and C3 are fermionic constants. Since \x 
is not arbitrary, this solution is not the same as that of £5. 



B Stationary solutions 

For the subalgebra £ 10 = {Pi + vQ 2 }, the obtained stationary solution is 

W = f{9 1 ) + m (0 2 -ux), P = g(9 1 ) + //A(0 2 - ux) (81) 

which contains two arbitrary functions, / and g, of two fermionic arguments, 9\ and 
9 2 — vx- This can be reduced to the form 

W = C + Ci6>i + m (9 2 -ux), P = C 2 + Cgfli + fiX(9 2 - ux) , (82) 
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where C and C 2 are bosonic constants and Cj_ and C3 are fermionic constants. This 
solution is linear in x, 9± and 9 2 . 

For the subalgebra £n = {Pi + eP 2 + V.Q2}, the obtained stationary solution is 

W = fie,) + ^9 2 -ux), P = s(0i) + A(# 2 - zat) (83) 

which contains two arbitrary functions, / and g, of one fermionic argument 9\. This can 
be reduced to the linear form 

W = C + C 1 9 1 + f±(9 2 -vx), P = C 2 + C 1 e 1 + \{e 2 -vx), (84) 

where Co and C 2 are bosonic constants and Cj_ and C3 are fermionic constants. 

For the subalgebra £ 13 = {P 2 + I^Qi}, the obtained stationary solution is 

W = f{6 2 ) + ^776*1 + fjjyx, P = g(6 2 ) + yXdx + [xkx (85) 

which contains two arbitrary functions, / and g, of one fermionic argument 9 2 . This can 
be reduced to the linear form 

W = C + Ci6 2 + fjjfii + ^vx, P = C 2 + Cs9 2 + fiXOx + \i_kx, (86) 

where Co and C 2 are bosonic constants and C\ and C3 are fermionic constants. Since /i 
is not arbitrary, this solution is not the same as that of £n. 

For the subalgebra £ 2 i = {Li + eP 2 }, the obtained stationary solution is 

W = x 3/2 (±ix~ 1/2 8i + inf(0 2 )), P = x^\ux- 1/2 9 1 + rMl{d 2 )) (87) 

which contains two arbitrary functions, / and g, of one fermionic argument 9 2 and is 
parametrized by the polytropic exponent 7. The solution can be expressed as 

W = ^x9 x + x z/2 r ]} C 1 9 2 + x 3/ \C , 

(88) 

P = x wf2 [y X -^e x + rj 2 C 1 9 2 + ryf 2 ) 
where C 2 is a bosonic constant and C\ and C 3 are fermionic constants. 

For the subalgebra £ 22 = \L\ + vQ 2 }, we obtain the stationary singular solution 

W = fj,x$i + x 3/2 ??# 2 - \x z/2 (In x)r\y_ 

( * \ ( 7+5 \ , 1 7+5 \ (89) 
P = (aA 2 1 +2)Q l _)_ x ^ 2 i+ 2 > p9 2 - ixv27+2J(lnx)pz/ 

which is parametrized by the polytropic exponent 7. 



C Solutions in the form of a reduced equation 

For the subalgebra £ 8 = {L 2 }, for the case where 7 = 3, we obtain the following reduced 
partial differential equations 

A + r 2 A T2 - 2A x A eiT2 + t 2 A x A 01 + 29 1 A X A XT2 - 9 1 t 2 {A x ) 2 = 0, (90) 
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where A is the superfield described in Table II with invariant r 2 = t 1 / 2 9 2 . The superfield 
W is expressed in terms of A through the relation 

W = r 1/2 A(x,e 1 ,T 2 ) (91) 

Here, P is a constant bosonic. This solution represents a non-stationary simple wave 
involving both bosonic and fermionic variables. 

For the subalgebra £i 2 = {L 2 + ePi}, for the case where 7 = 3, we obtain the following 
symmetry reduction 



-A- eA c - r 2 A T2 + 2A a A dlT2 - r 2 A a A ei - £T 2 A a A a9l - 29 1 A a A, 
+ Oir 2 (A a ) 2 + eB^AaAaa = 0, 



where A is the superfield described in Table II with invariants a = x — |e In t and r 2 = 
t~ x l 2 Q 2 . The superfield W is expressed in terms of A through the relation 

W = r 1 / 2 A(a,9 1 ,T 2 ) (93) 

Here, P is a constant bosonic. This solution represents a non-stationary simple wave 
involving both bosonic and fermionic variables. 

For the subalgebra £i 6 = {L 2 + /iQi}, for the case where 7 = 3, we obtain the following 
symmetry reduction 

- A + flTiAcr - flA Tl - T 2 A T2 + 2A a A TlT2 - T 2 A a A Tl + HT 2 (A a ) 2 + llTx^AaAa^ / 

- — — (94) 

- 2T 1 A a A aT 2 + 7iT 2 (At) 2 + Vn^AaAan =0 

where ^4. is the superfield described in Table II with invariants a = x + |/i#iln£, T\ = 
9 1 — iyulnt and r 2 = t~ 1 l 2 6 2 . The superfield W is expressed in terms of A through the 
relation 

W = r 1/2 A(a,T l ,r 2 ) (95) 

Here, P is a constant bosonic. This solution represents a non-stationary simple wave 
involving both bosonic and fermionic variables. 

For the subalgebra £ 20 = {L 2 + sP 1 + fiQi}, for the case where 7 = 3, we obtain the 
following symmetry reduction 

- A + \rr\Av - eAo - [jA t1 - t 2 A T2 + 2A a A T1T2 - r 2 A a A T1 + ^T 2 {A a ) 2 

+ fiTi t 2 A C A C T\ - er 2 A a A aTl - 2T 1 A a A aT2 + nT 2 (A a ) 2 + et^A^A^ (96) 

+ \vr\T 2 A c A CTl = 

where ^4 is the superfield described in Table II with invariants a = x + lnt — \e lnt, 

ri = ^1 — |/xlnt and r 2 = t~ l l 2 Q 2 . The superfield is expressed in terms of A through 
the relation 

W = r 1/2 A(a,r 1 ,r 2 ) (97) 

Here, P is a constant bosonic. This solution represents a non-stationary simple wave 
involving both bosonic and fermionic variables. 
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For the subalgebra £ 24 = {L\ + kL 2 }, for the case where 7 = 2 and k = 7, we obtain 
the following symmetry reduction 

-4A- 2oA a - nA T1 - 7t 2 A T2 + UA T1T2 A a - Ar 2 A T1 A a - 2ot 2 A utx A - r 2 A T1 A a 

- UnAa-^Aa + 4r 1 T 2 (A) 2 + 2<JT 1 T 2 A aa A a = 

(98) 

_i j_ 

where A is the superfield described in Table II with invariants a — t 7x, T\ — t and 
_i 

r 2 = t 2 2 . The superfield W is expressed in terms of A through the relation 

W = r$A(a,T 1 ,T 2 ) (99) 

Here, P is a constant bosonic. This solution represents a non-stationary simple wave 
involving both bosonic and fermionic variables. 

D Wave solutions 

For the subalgebra £ 3 = {Pi + nQ±}, the reduced system of equations is 

P t - 2/xP T1 P Tie2 - 2/xP T1 fl 2 P tT1 + P Tld2 W x + 9 2 P tT1 W x + mP Tld2 W x + vr x e 2 P tTX W x = 0, 
W t + 2W Tld2 W x + 26 2 W tT1 W x + I^W^W* + 2^Q 2 W tT ,W x = 

(100) 

and the obtained solution is 

W = [ itf 1 (r 1 ,r 2 ) + f 2 ( Tl ,T 2 ), P = fd 9l (r U T 2 )+g 2 (T U T 2 ), (101) 

where fi and g 1 are fermionic-valued functions and f 2 and g 2 are bosonic-valued functions 
of the two arguments 

T\—0\ — fix and r 2 = 9 2 + /it. (102) 
The functions f±, f 2 , g\ and # 2 satisfy the relations 

fi — 2 /2,-n 7^/2,1-1 — /2,72) #1 = 92,TiT 2 f2,Ti + 2g 2yT1 g 2)T1T2 — g 2tT2 , (103) 

where / 2 and g 2 are completely arbitrary functions of t\ and r 2 . This solution has a 
degree of freedom of two arbitrary functions of two variables which combine the fermionic 
variables 6\ and 9 2 with the bosonic variables x and t. It represents a propagation wave 
in both the bosonic and fermionic variables. Since the rank of the matrix 

~ dW dW dW dW 

dx dt 80! 002 

dP dP dP dP 

_ dx dt d0! de 2 

is 2, this solution is called a Riemann scattering double wave. For the specific case where 
we choose 

j2 = e kT1T \ 92 = e lT ^, (105) 
we obtain the propagating wave solution 

W = fit (2k 2 9 2 + kOi) + 1 + kO^t + k6 2 fix + kQ x d 2 , 

(106) 

p = fit (ki + 21 2 ) e 2 + (die! + 1 + w^t + w 2 fix + iq x q 2 



(104) 
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For the subalgebra £4 = {Li}, we have obtained the explicit stationary solution in 
monomial form 

W = x 3/2 fi6 2 , P = Jj$*)v0 2 (107) 
which depends on the polytropic exponent 7. 

For the subalgebra £7 = {P2+KQ2}, the obtained non-stationary simple wave solution 
involving both bosonic and fermionic variables is 

W = g{fd + k6 1 + f±ue 2 ), P = f(jrt + k6 1 + ±v6 2 ) (108) 

which involves a bosonic constant k and two arbitrary functions, / and g, of one fermionic 
argument [Mt + k9\ + \ivQ 2 . 

For the subalgebra £ 23 = {L± + eP 2 + vQ 2 }, the obtained solution is the simple wave 

W = ^ ,2 A(a) P = rp(£&)a (109) 

which contains one arbitrary function A of one argument a = t + \vQ 2 In x — In x and 
is parametrized by the polytropic exponent 7. In the specific case where one chooses 
A(a) = e~ ka , we obtain 

W = fix^e~ kt (l + In x) , P = rpc^a. (110) 
Here, W involves damping in time t. 

The remaining subalgebras have invariants which possess a non-standard structure in 
the sense that one does not obtain standard symmetry reductions from them (as described 
irP^). Such subalgebras are distinguished by the fact that each of them admits an invari- 
ant expressed in terms of an arbitrary function of the superspace variables, multiplied by 
a fermionic constant. We list in Table III the invariants expressed in terms of an arbitrary 
function of the superspace variables for each case. 



Table III: List of subalgebras with nonstandard invariants. In each case, / is an arbitrary 
function of its arguments. 



Subalgebra 


Non-standard invariant 


£2 = {/uQi} 




£ 6 = {V.Q2} 


vj (x,t,e 1 ,0 2 ,w,p) 


£14 = {fj-Qi +EQ2} 


iivj (x,t,9i,e 2 ,W,P) 


£15 = {^2+^1+^2} 


iivj (x,e u e 2 ,w,p) 


£is = {-Pi + mQi + EQ2} 


fj, E f (t,9 u e 2 ,w,p) 


£19 = {Pi + eP 2 + vQi + EQ2} 


HvJ(9 x ,0 2 ,W,P) 



V Concluding remarks and future outlook 

The objective of this paper was to construct a supersymmetric extension of the polytropic 
gas dynamics equations through a superspace involving two independent variables and two 



23 



bosonic superfields and to perform a systematic group-theoretical analysis of this exten- 
sion. We found a Lie superalgebra of symmetries of the proposed supersymmetric model 
pip which included two translations, two dilations and two supersymmetric transforma- 
tions. Decomposing this superalgebra into a direct sum, we used the Goursat method 
to classify its one-dimensional subalgebras into 24 conjugacy classes under the action of 
the associated supergroup. By using the symmetry reduction method, we determined the 
invariant solutions of the supersymmetric model under consideration. By postulating a 
wave-like form for the solutions, we obtained a number of classes of solutions involving 
arbitrary functions of one or two variables. In the particular case where the arbitrary 
functions involved only 9\ and 82, we were able to explicitly find the most general form 
of the solution through a truncated Taylor expansion. We constructed algebraic-type so- 
lutions and propagation waves (including simple and double waves) in explicit form. We 
conclude that the solutions of the supersymmetric polytropic gas dynamics model (TTTT) 
include propagation waves and their superpositions involving both bosonic and fermionic 
variables. This phenomenon is analogous to that of the multisolitonic wave solutions 
obtained in the case of supersymmetric extensions of integrable models. 

The following question can be considered. The freedom of the obtained solution can 
be found from the Cauchy data when t = 0. However, the aspect of the Cauchy problem 
involving uniqueness and continuous dependence on initial conditions remains an open 
problem for the case of hydrodynamic-type equations. This will be addressed in a future 
work. 
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